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ABSTRACT 

We present a full-sky model of polarized Galactic microwave emission based on three years of 
observations by the Wilkinson Microwave Anisotropy Probe (WMAP) at frequencies from 23 to 94 
GHz. The model compares maps of the Stokes Q and U components from each of the 5 WMAP 
frequency bands in order to separate synchrotron from dust emission, taking into account the spatial 
and frequency dependence of the synchrotron and dust components. This simple two-component 
model of the interstellar medium accounts for at least 97% of the polarized emission in the WMAP 
maps of the microwave sky. Synchrotron emission dominates the polarized foregrounds at frequencies 
below 50 GHz, and is comparable to the dust contribution at 65 GHz. The spectral index of the 
synchrotron component, derived solely from polarization data, is -3.2 averaged over the full sky, with 
a modestly flatter index on the Galactic plane. The synchrotron emission has mean polarization 
fraction 2-4% in the Galactic plane and rising to over 20% at high latitude, with prominent features 
such as the North Galactic Spur more polarized than the diffuse component. Thermal dust emission 
has polarization fraction 1% near the Galactic center, rising to 6% at the anti-center. Diffuse emission 
from high-latitude dust is also polarized with mean fractional polarization 0.036 ± 0.011. 

Subject headings: polarization, cosmic microwave background, radio continuum: ISM, dust 



1. INTRODUCTION 

The Wilkinson Microwave Anisotropy Probe {WMAP) 
has mapped the full sky in the Stokes /, Q, and U pa- 
rameters on angular scales 6* > 0?2 in 5 frequency bands 
centered at 23, 33, 41, 61, and 94 GHz, denoted K , 
Ka, Q, V, and W respectively (jBennett et al.l l2003bl : 
iHinshaw et al.l 120071 ). Polarized emission at microwave 
frequencies is dominated by foreground emission from 
the Galaxy, which at degree angular scales is brighter 
than the cosmolo gical signal in all of the WMAP bands 
(jPage et al.l 1200 7^. Extraction of cosmological informa- 
tion from the WMAP polarization data requires fitting 
for the astrophysical foregrounds. 

Foreground polarization at microwave frequencies is 
dominated by a superposition of synchrotron and ther- 
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mal dust emission. Synchrotron emission results from 
the acceleration of cosmic ray electrons in the Galactic 
magnetic field. For a power-law distribution of electron 
energies N(E) cx E~p propagating in a uniform mag- 
netic field, the resulting emission is partially polarized 
with fractional linear polarization 



p+1 
p + 7/3 



(1) 



aligned perpendicu lar to the magnetic field 

(jRvbicki fc LightmanI I1979D . The frequency depen- 
dence of synchrotron emission is also related to the 
electron energy distribution, T{i/) cx ly^ with spectral 
index 

/3 = -^ (2) 

where T is in units of antenna temperature. For spec- 
tral index /? « —3 observed at microwave frequencies, 
synchrotron emission could have fractional polarization 
as high as fs ^ 0.75. Line-of-sight and beam averag- 
ing effects will tend to reduce the observed polarization 
by averaging over regions with different electron energy 
distribution or magnetic field orientation. 

Thermal dust emission can be partially polarized as 
prolate dust grains align wi th their long axis perpen - 
dicular to the magnetic field (|Davis fc Greensteinl[l95lh . 
The dust emission efficiency is greatest along the long 
axis, leading to partial linear polarization perpendicular 
to the magnetic field (thus following the same direction 
as synchrotron emission) . The fractional polarization de- 
pends on the grain size distribution and is typically a 
few percent at millimeter waveleng ths (Hildcbrand 198^; 
'Hilde brand et al.lll999l : IVaillancourtll2002l . 

Other sources of polarized emission are possible. Ex- 
tragalactic radio point sources, although detected in 
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Fig. 1. — Estimated signal variance at each WMAP frequency 
band for 3? 7 pixels. The solid line shows a two-component fit with 
power-law synchrotron and thermal dust emission assuming no spa- 
tial correlations between the two components. Allowing spatial 
correlations between the synchrotron and dust better reproduces 
the observed spectra (dotted lines). The thick dashed line shows 
the 68% confidence noise level. Note the change in scale between 
panels: the top panel shows high-latitude data outside the P06 sky 
mask, while the bottom panel shows data inside P06 dominated by 
emission from the Galactic plane. 



the unpolarized temperature maps, contribute neghgi- 
bly to the observ ed polarization on degree angular scales 
(jPage et al.ll2007D . Non-thermal emis sion from dust via a 
popu lation of small spinning grains ( Draine fc LazarianI 
1998) or partially magnetized grains ( Draine fc LazarianI 
1999f ) has been proposed to explain the observed correla- 
tion between far-infrared and microwave emission in the 
interstellar medium. While emission from spinning dust 
is thought to be largely unpolarized, emission from par- 
tially magnetized grains could be substantially polarized 

rovide d the grains consist of a single magnetic domain 

Draine & Lazarian 1999). 

iPage et al.l (|2007D present a simple model of foreground 
polarization based on the three-year WMAP data. Spa- 
tial templates for the synchrotron and dust emission, fit 
to the high-latitude sky, reduce the residual foregrounds 
in the corrected maps to levels below the polarization 
of the cosmic microwave background. This paper com- 
pares the WMAP polarization data to the unpolarized 
synchrotron and dust emission to model the polarized 
foregrounds, inchiding both the spatial and frequency 
variation of the foregrounds over the full sky. 



2. FREQUENCY SPECTRUM 

The signal variance at each WMAP frequency band 
provides a simple, model-independent estimate of the 
frequency spectrum of polarized emission. The data is 
a superposition of signal (either cosmic or foreground) 



(3) 



where i is a pixel index, a and b are map indices, a de- 
notes the Stokes Q oi U maps, and w is a pixel weight. 
We use different weights in different regions of the sky. 
Pixels at low Galactic latitude are signal-dominated and 
use unit weight, Wi = 1. At high latitudes, instru- 
ment noise becomes appreciable and we use noise weight, 
Wi = Ni/oQ where Ni is the number of observations in 
pixel i and ctq is the noise per observation. If the vari- 
ance is computed from a single map (a = b), instrument 
noise contributes a positive contribution to each pixel, 
(T"^) — ((s + n)^) = + {ri^)- We avoid this noise bias 
by computing the cross variance using maps from differ- 
ent years, a ^ h. Uncertainties in the noise model then 
affect only the uncertainty in the estimated variance but 
not the mean value itself. We estimate the uncertainty in 
the computed cross-variance using Monte Carlo simula- 
tions consisting of a model signal phis noise realizations 
derived from a Cho lesky decomposition of the noise ma- 
trices for each map parosik et al.|[2007l) . The noise simu- 
lations thus include the contribution of 1// noise as well 
as correlations between different pixels and between the 
Stokes Q and U maps. We add the same simulated fore- 
ground to each noise realization prior to computing the 
variance, to account for the additional scatter from one 
realization to the next caused by the signal-noise cross 
term 2sn. 

Figure [1] shows the frequency spectrum of the signal 
variance for two regions of the sky defined by the P06 
polarization mask (Pag e et al.ll2007l ): the "high signal" 
region inside the P06 mask, dominated by the Galactic 
plane and North Galactic spur where signals are clearly 
visible in the maps, and a "low signal" region outside 
the P06 mask where noise is more prominent. Both re- 
gions are dominated by a power-law component T (x 
with spectral index j3s ~ —3 consistent with synchrotron 
emission. The rise in variance from V band to W band 
(61 to 94 GHz) requires a second component with spec- 
tral index /^^ « 2 consistent with thermal dust. The 
observed variance is brighter than the expected cosmic 
signal (s^ « 0.14 /xK^ within 3? 7 pixels) at all frequen- 
cies. Detection of a cos mic signal requires subtraction of 
foreground polarization (jPage et al.ll2007l ). 

The observed steepening of the cosmic ray energy spec- 
trum is expected to steepen the synchrotro n emission 
spectrum relative to a pure power law (Miiller fc Tand 
fl987HV5ia[l989l) . The polarization variance does not 
show the expected steepening above 30 GHz, but rather 
shows a spectrum flatter than a pure power law. This 
could indicate either the presence of a third emission 
component or non-zero spatial correlations between the 
synchrotron and dust components. If the synchrotron 
and dust components each follow a power law in fre- 



quency, Tsynch 
ance becomes 

= S' 



Sv^ and T, 



dust 



Dv^ , then the vari- 



2/3d 



4- 2rSDu'^^+'^-' 



(4) 



where r is the spatial correlation between the two com- 
ponents. We fit the observed variance data to a 2- 
component model including spatial correlation, and find 
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a broad minimum with r ~ 0.8 (dotted lines in Fig- 
ure [T]). Since the unpolarized synchro tron and dust 
foregr o unds tend to b e highly correlated (|Bennett et al.l 
l2003al: iHinshaw et al.l [20071. it is not surprising to find 
that the correlation persists in polarized emission from 
the same media. 

3. FOREGROUND POLARIZATION 

The WMAP 3 year data map the polarized synchrotron 
emission at high signal to noise ratio. Polarized dust 
emission is less well constrained, particularly at high 
Galactic latitudes. We may use maps of the unpolarized 
dust to constrain the maximum amplitude of polarized 
dust emission. The dust polarization is 

Pd{n, v) = fd{n)Td{h, y) (5) 

where P = (Q^ + f/^)"-^ is the polarized amplitude, Td 
is the unpolarized amplitude, and fd is the fractional 
polarization in direction fi. Since fd must lie in the range 
[0,1], the polarized dust amplitude can not exceed the 
total unpolarized emission in any pixel. 

We thus model the polarized foregrounds as a super- 
position of synchrotron and thermal dust emission, 

Q{h, iy)^Ps{n)S{n, v) cos(27s(n)) 

+ fd{n)Td{n)D{v) cos(27d(n)) 

U{h, v) = P, {n)S{n, v) sin(27, {h)) 

+ fd{h)Td{h)D{v)^iu{2-fd{h)) (6) 

where Ps{n) is the (polarized) synchrotron amplitude. 
We model the frequency dependence of the dust as a 
spatially invariant power law in antenna temperature, 

Di.)={^f (7) 

with spectral index fid = Synchrotron emission may 
be approximated as a power law, but the spectral index 
varies with position on the sky and may steepen with 
frequency. We thus model the synchrotron frequency de- 
pendence as 

/ \(}s{n)+Clog.{v/uo) 

5(^,n)=(^-j (8) 

allowing the spectral index I3s{n) to vary with position. 
The model includes a logarithmic steepening of the spec- 
trum with frequency, with C ~ —0.4 typical of the syn- 
chrotron e mission in the WMAP unpolarized tempe ra- 
ture data (|Hinshaw et al.ll2007t [Bennett et al.ll2003aD . 

The polarization angle 7(n) = 0.5 arctan([/(n)/(5(n)) 
is defined with respect to the Galactic meridian. Syn- 
chrotron emission dominates the lowest frequency chan- 
nels, providing a template for the synchrotron polariza- 
tion angle jsin). We use the 22 GHz (K-band) data 
alone to compute the synchrotron polarization angle 7^. 
In principle, this contains a small level of contamina- 
tion from the cosmic signal. An alternative method com- 
putes the polarization angle from the K-Ka or K-Q dif- 
ference maps. This eliminates any cosmic signal at the 
cost of increased noise. We use Monte Carlo simulations 
to compare the two methods. For noise levels typical of 
the WMAP three-year polarization maps, direct fitting 
to the K-band data better reproduces the input pattern 
used to generate the simulations. 




Fig. 2. — Spectral index of polarized synchrotron emission be- 
tween 22 and 33 GHz (MoUweide projection in Galactic coordi- 
nates). The map has been smoothed with a position-dependent 
tophat varying from 7° radius near the plane to 18° radius at high 
latitude. 

In this paper, we do not assume that the dust follows 
the same polarization direction as the synchrotron, so 
that the synchrotron polarization angle 7^ (n) may differ 
from the dust polarization angle jdin). This allows for 
dust emission to originate from a different spatial dis- 
tribution compared to synchrotron emission along each 
line of sight. We adopt the dust polarization angles from 
the starhght model described in lPage et al.l (|2007[ ). and 
compare this model to one in which 7^ = 7., everywhere 
on the sky. 

The foreground model does not assume that the syn- 
chrotron spectral index is identical along all lines of 
sight, but fits a map of the spectral index as described 
in ^ 13. II Given the polarization angles 7s (n) and ^din) 
and the map of the synchrotron spectral index /3s (n) , we 
fit the WMAP 3- year band-averaged polarization maps 
to derive the polarized synchrotron amplitude Ps(n), 
synchrotron spectral correction C, and the dust frac- 
tional polarization fdin) (Eq. [S]). Based on signal to 
noise, we fit these parameters simultaneously but on dif- 
ferent angular scales. Synchrotron emission dominates 
the low-frequency data. We fit for the synchrotron am- 
plitude Ps{n) (normalized to K band) in each of 3072 
equal area pixels at HEALPix resolution 4 or A^side = 16 
(jGorski et al.ll2005[ ). Polarized dust emission is consid- 
erably fainter than synchrotron in the WMAP frequency 
bands. We thus fit for the dust fractional polariza- 
tion fdin) in 6 coarse sky regions, holding fd constant 
within each region while simultaneously fitting for the 
synchrotron amplitude in each of the smaller HEALPix 
pixels. We fit for the synchrotron spectral correction 
C in two regions defined by the P06 polarization mask: 
one value for the high-latitude region outside the mask 
and an independent value for the low-latitude region in- 
side the mask. W e model the unpolarized dust emission 
Td{n,v) using the iFinkbeiner et all (|1999D (FDS) model 
8, evaluated at 94 GHz and scaled with uniform spectral 
index j3d = The results do not depend sensitively on 
the choice for Pd- 

3.1. Polarized Synchrotron Model 
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Fig. 3. — Spectral index of polarized synchrotron emission binned 
by Galactic latitude. The spectral index steepens off the Galac- 
tic plane. Similar behavior is observed in external edge-on spiral 
galaxies. 

Synchrotron emission may be approximated as a power 
law in antenna temperature, with spect ral index Psjn) 
varyin g with position across the sky. iHinshaw et al.l 
([2007") present a maximum-entropy model of unpolar- 
ized synchrotron emission derived from WMAP data. We 
may use this maximum-entropy model to specify the fre- 
quency dependence of the polarized synchrotron emis- 
sion, including both the spatial variation and the steep- 
ening of the spectrum with respect to a pure power law. 
The unpolarized maximum-entropy model, however, suf- 
fers from confusion with free-free emission, particularly 
at low Galactic latitudes. 

Synchrotron emission dominates the polarization data 
at 22 and 33 GHz (K and Ka bands), with negligible con- 
tribution from thermal dust (Fig. [J). We thus use the 
WMAP polarization data at K and Ka bands to map the 
synchrotron spectral index j3s{n). We fit for the spectral 
index in each of 3072 pixels, imposing a uniform prior 
—4.0 < Ps < —2.0 to reduce the effects of instrument 
noise. The resulting map is still noise dominated, partic- 
ularly at high latitude. We further reduce the noise con- 
tribution by smoothing the spectral index map using a 
tophat window function. We select the tophat radius and 
verify the accuracy of the fitting procedure using Monte 
Carlo realizations consisting of a fixed Galactic signal 
added to random realizations of instrument noise. At 
each pixel, we smooth the spectral index maps from each 
realization using a progressively larger tophat radius, un- 
til the standard deviation of the ensemble of smoothed 
realizations, evaluated at that pixel, falls below a thresh- 
old 6f3s{n) < 0.1. The resulting smoothing radius varies 
across the sky from 7° near the plane to 18° at high 
latitude. 

We investigate the dependence of the recovered spec- 
tral index map on the prior imposed during the pixel-by- 
pixel fit. Changing the prior to a wider range —4.5 < 
Ps < -1.5 or a narrower range —3.5 < Ps < —2.5 pro- 
duces a systematic shift 5 (3s ~ 0.1 in the mean of the 
smoothed spectral index map, comparable to random un- 
certainty, but does not significantly alter the structure on 
smaller angular scales. Our choice of prior is informed by 
Monte Carlo realizations, where the ensemble average of 
the smoothed realizations accurately recovers the noise- 



less spectral index map used to generate the simulations. 

Figure [2] shows the polarized synchrotron spectral in- 
dex from the WMAP 3- year data. The mean value aver- 
aged over the full sky is {(3s) = —3.2. The spectral index 
is flatter along the Galactic plane with steeper values at 
higher latitude. Figure [3] shows the polarization spectral 
index binned by Galactic latitude. 

The smoothing applied to the spectral index maps cre- 
ates spatial correlations in the smoothed maps. We ac- 
count for this in Figure [3] by plotting the effective un- 
certainty in the binned data, SP^^ = at/y/ N°^, where 
(Tj is the standard deviation of the smoothed values from 
all pixels in the i^^ latitude bin, and Nf^ is the num- 
ber of independent spatial regions within that bin. The 
spectral index steepens from Ps ~ —3.05 along the plane 
to Ps « —3.25 at latitude |6| > 30°. Polarization data 
provide a measurement of the synchrotron spectral index 
independent of the temperature data and not subject to 
confusion from competing emission mechanisms. 

Figure [4] shows the model polarized synchrotron emis- 
sion. The polarized synchrotron amplitude Ps (n) is dom- 
inated by the lowest frequency data and is nearly iden- 
tical to the K-band map. We derive the synchrotron 
fractional polarization by comparing the K-band model 
amplitude to the unpolarized K-band maximum-entropy 
synchrotron model (,Hinshaw et al.,i2007. ). assuming that 
this model accurately represents the unpolarized syn- 
chrotron amplitude. Although the Galactic plane domi- 
nates the polarized emission in the maps, the correspond- 
ing fractional polarization is low, indicative of contribu- 
tions from multiple emitting regions along the line of 
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Fig. 4. — Polarized synchrotron model parameters. (Top) An- 
tenna temperature P = (Q^-t-f/^)"-^ of polarized synchrotron emis- 
sion at K band. Vectors indicate the polarization direction. (Bot- 
tom) Synchrotron fractional polarization, derived by dividing the 
polarized model by the unpolarized maximum-entropy synchrotron 
model at K band. The polarization fraction is 3-5% on the Galac- 
tic plane, increasing above 20% within the North Galactic Spur 
and its southern extension. 
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Fig. 5. — Histogram of synchrotron fractional polarization fs- 
(Top) High latitude region outside the P06 mask. (Middle) Mid- 
latitude region inside P06 but with |b| > 5°. (Bottom) Galactic 
plane, |fe| < 5°. 

sight. The North Galactic Spur and its southern exten- 
sion are more prominent in polarization than in inten- 
sity, indicating localized enhancements in the polariza- 
tion fraction. Figure [5] shows a histogram of the syn- 
chrotron fractional polarization for several regions of the 
sky. The Galactic plane < 5°) is less than 10% po- 
larized, with mode fs — 0.02 and mean fs — 0.05. The 
North Galactic Spur has enhanced fractional polariza- 
tion, fs ~ 0.3, while the high-latitude sky outside the 
P06 mask has mode fs — 0.14 and mean fs — 0.15. 

Figure [6] shows the fractional polarization binned by 
Galactic latitude. Within narrow latitude bins, the dis- 
tribution of observed fractional polarization values ap- 
proximates a Gaussian distribution. Figure [6] shows the 
mean (fi) within the i^^ bin with uncertainty Sfi — 
o'i/V^i defined by the standard deviation ai of the Ni 
pixels within that bin. The uncertainties thus reflect 
the precision to which the mean may be defined within 
each latitude bin, which average over a larger range of 
fractional polarization (c/ Figure [5]). The mean frac- 
tional polarization rises linearly from 4% near the plane 
to 20% at latitude \b\ > 50°. The north polar cap region 
is more polarized than the south polar cap, although 
this may be an artifact of the unpolarized synchrotron 
model used to define the synchrotron fractional polar- 
ization. Both the low polarization in the plane of the 
Galaxy and the trend toward higher fractional polariza- 
tion at high latitude are similar to previous measure- 
ments of external spiral galaxies. Edge-on spirals typ- 
ically exhibit 1-3% polarization in the disk, rising lin- 
early to 10 -20% at dis tances of several kpc from the 
ma.ior axis jjHummel efa l. 1991; Sukumar & Allen 1991; 
Dumke et al.lll995HDui^ke fc Krauselll998HDumke et all 
2000() . The large-scale polarization of the Milky Way is 
thus unremarkable compared to other galaxies. 



The foreground model directly fits the polarized syn- 
chrotron amplitude Ps (n) and compares this fitted ampli- 
tude to a model of the unpolarized synchrotron emission 
to estimate the fractional polarization fs . A recent com- 
parison of the unpolarized WMAP K-band map at 22 
GHz with a full-sky map at 19 GHz suggests that much 
of the emission near the Galactic plane results from a 
component spatially correlated with thermal dust emis- 
sion but with with spectral index consistent 
with spinning dust (|Boughn fc Pobei1l2006D . To the ex- 
tent that the maximum-entropy synchrotron model over- 
estimates the unpolarized synchrotron amplitude by ig- 
noring spinning dust, the fractional polarization fs de- 
rived above becomes a lower limit to the true synchrotron 
fractional polarization. Since the values for fs derived 
using the maximum-entropy model agree with observa- 
tions of external edge-on spiral galaxies at much lower 
frequencies (where spinning dust is negligible), requiring 
a substantial fraction of the unpolarized emission at 22 
GHz to originate from spinning dust would require some 
compensating depolarization of the synchrotron compo- 
nent {e.g. line-of-sight effects or a tangled magnetic field) 
to keep the Milky Way in family with similar external 
galaxies. 

3.2. Polarized Dust Model 

Polarized dust emission is brighter than synchrotron 
only in the highest frequency channel, and only signifi- 
cantly so on the Galactic plane. Since the dust fractional 
polarization fd is by definition positive, noise in the W 
band data can create a positive bias in the fitted dust 
solution. We estimate this bias using Monte Carlo sim- 
ulations, comparing the recovered synchrotron and dust 
parameters to known inputs in each coarse sky region for 
300 realizations consisting of synchrotron, dust, CMB, 
and instrument noise generated using Ch qlesky decom- 
posit ion of the full noise matrix in each DA (jJarosik et al.l 
l2007f) . The recovered dust solution is nearly identical to 
the input, with bias less than 10% of the statistical un- 
certainty. We do not apply a correction for this effect. 

Table [T] shows the fractional polarization of the dust 
in each coarse sky region. The fractional polarization in 
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Fig. 6. — Synchrotron fractional polarization fa binned by Galac- 
tic latitude for the northern (filled symbols) and southern (open 
symbols) hemispheres. Each point shows the mean fractional po- 
larization fs for all pixels in a latitude bin with signal-to-noise ratio 
greater than one. The error bars show the uncertainty in the mean 
based on the scatter of values within each latitude bin. 
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TABLE 1 

Dust Fractional Polarization 

Region Galactic Galactic Fractional 



Latitude Longitude Polarization 

(dcg) (dcg) /d 

T -10 to 10 -135 to -45 0.010 ± 0.003 

2 -10 to 10 -45 to 45 0.010 ± 0.005 

3 -10 to 10 45 to 135 0.012 ± 0.005 

4 -10 to 10 135 to -135 0.064 ± 0.013 

5 |b| > 10, inside P06 0.020 ± 0.014 

6 Outside P06 0.036 ± 0.011 



the plane increases away from the Galactic center, with 
fd = 0.01 at the Galactic center rising to 0.06 at the 
anti-center. The fractional polarization increases with 
latitude, reaching a value fd — 0.036 ± 0.011 outside the 
P06 polarization mask. 

Equation [6] fits a single parameter to the high- latitude 
dust (outside the P06 mask). We have repeated the fit 
using more "coarse" sky regions outside the P06 mask. 
The faint amplitude of the dust signal combined with the 
lack of small-scale structure in the high-latitude template 
map D{h) prevent a statistically significant determina- 
tion of substructure in the high-latitude polarized dust 
signal. 

4. DISCUSSION 

We model the polarized foreground emission in each 
pixel as a superposition of synchrotron emission with 
spatially variable spectral index determined from the 
WMAP polarization data plus dust emission traced by 
FDS dust model scaled by a spatially variable fractional 
polarization. Figure [7] shows the model and residuals for 
each frequency band. A simple test compares the vari- 
ance of the foreground-cleaned maps to the variance of 
the uncleaned maps (Equation[3]). The foreground model 
removes over 97% of the power at Ka and Q bands where 
foregrounds are brightest. Outside the P06 mask, the 
variance of the cleaned maps is close to the value ex- 
pected for instrument noise alone. 

Table [H c ompa res the template model described in 
iPage et aD ()2007l ) to the pixel model described in this 
paper, using goodness-of-fit statistic 

= T,ijal3RiaN~.lpRji3 (9) 

where N^^ is the inverse noise matrix and Ria is the 
(data- model) residual in pixel i and Stokes parameter a. 
The synchrotron amplitude Ps and polarization angle 7s 
in the pixel model are derived almost entirely from the 
K-band data, leaving nearly zero residuals at K band. 
Since the synchrotron amplitude and polarization angle 
may be treated as an external (K-band) template, we as- 
sess the goodness of fit using only the 4 higher-frequency 
bands, counting only the parameters fitted using data at 
these bands. The high-latitude region outside the P06 
mask contains 2267 pixels in each of the Stokes Q and U 
maps in each of the 4 frequency bands. The synchrotron 
spectral index map has 36 independent regions outside 
the P06 mask (after accounting for smoothing), requir- 
ing 36 parameters in the extended model. We add one 
additional parameter for the synchrotron spectral correc- 
tion C plus another for the dust fractional polarization 
fd to obtain a total of 38 parameters fitted to 18136 data 
points outside the P06 mask. For the 805 pixels inside 



the P06 mask we fit 31 parameters for the synchrotron 
spectral index, 5 parameters for dust, plus one param- 
eter for the synchrotron spectral correction for a total 
of 36 parameters fitted to 6440 low-latitude data points. 
The template model, by comparison, fits only 2 parame- 
ters to each band, one each for the dust and synchrotron 
template maps. 

Fitting a spatially variable synchrotron spectral index 
provides a statistically significant improvement (Ax^ = 
447 for 30 additional model parameters outside the P06 
mask) compared to the template model, which assumes 
the same spatial pattern for the synchrotron at all bands. 
Most of this improvement comes at Ka and Q bands 
where synchrotron is brightest. Although the pixel-based 
model has a lower a-t V band, the difference at this 
band alone is not statistically significant. The template 
model has slightly better at W band outside the P06 
mask. Not surprisingly, the pixel-based model provides 
a better fit at all bands inside the P06 mask. 

The pixel-based model finds a synchrotron spectrum 
marginally flatter than a power law, with best-fit spec- 
tral correction C = 0.3 ±0.1. This result is unchanged 
with modifications of the model (polarization angle, syn- 
chrotron spectral index map, dust spectral index) and 
is independent of the sky cut. The apparent flattening 
is a robust feature of the data, derived independently 
from analyses of the data variance (jj 2), the pixel-base d 
model (§3), and the template model (|Page et al.l 120071 ). 
We obtain the same value for C fltted independently to 
the regions inside and outside the P06 mask. We note 
that the value for C is heavily weighted by the data at Q 
and V bands where both the synchrotron and dust are 
faint. Forcing C = to fit a pure power law increases 
the by 10 outside the P06 mask, while a steepening 
spectrum (C = —0.4) increases by 45. In both cases 
the the fitted dust parameters shift slightly compared to 
the best-fit model. We define a systematic uncertainty 
for the dust fractional polarization as the change in fd 
induced as the synchrotron spectral correction parame- 
ter C varies from 0.3 to -0.4. The uncertainties reported 
in Table [T] include this systematic uncertainty added in 
quadrature with the statistical uncertainty. 



TABLE 2 
OF Residual Maps 



Model 


Band 


r 


DOF 


X^/DOF 


Pixel Model 


Ka 


4604 


4524 


1.018 


Outside P06 


Q 


4521 


4524 


1.000 




V 


4600 


4524 


1.017 




W 


4774 


4524 


1.055 




All 


18505 


18098 


1.023 


Pixel Model 


Ka 


1626 


1601 


1.016 


Inside POO 


Q 


1541 


1601 


0.963 




V 


1670 


1601 


1.043 




W 


1773 


1601 


1.107 




All 


6611 


6404 


1.032 


Template Model 


Ka 


4968 


4532 


1.096 


Outside P06 


Q 


4616 


4532 


1.019 




V 


4608 


4532 


1.017 




w 


4760 


4532 


1.050 




All 


18952 


18128 


1.045 


Template Model 


Ka 


2403 


1608 


1.494 


Inside P06 


Q 


2081 


1608 


1.294 




V 


2704 


1608 


1.682 




W 


4946 


1608 


3.076 




All 


12034 


6432 


1.871 
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Fig. 7. — Polarized foreground model and residuals for the pixel-by-pixel polarization model. (Left) Antenna temperature of polarized 
emission at each band. The polarization direction is not shown. (Right) Residuals after subtracting foreground model. For visualization, 
the residuals P = \/Q2 + 1/2 

are shown, corrected for the positive noise bias using Monte Carlo simulations. The model and x 
computed directly from the Stokes Q and U maps. 
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The analysis quantifies the dependence of the syn- 
chrotron spectral index map on the choice of prior used 
in the fitting process (§3.1). We repeat the full fore- 
ground analysis, replacing the spectral index map I3s{n) 
in Eqs. [6] - [8] with a new map derived using a differ- 
ent prior before re-computing the best-fit values for the 
foreground model. We find a broad minimum in near 
prior —4.0 < (3s < —2.0, consistent with the Monte Carlo 
tests described in §3.1. Changing to a more restrictive 
prior —3.5 < (3s < —2.5 increases by 33, while the 
looser prior —4.5 < f}s < —1.5 increases x^ by 2. 

A x^ test confirms that the structure in the syn- 
chrotron spectral index map is dominated by a steep- 
ening of the spectral index off the Galactic plane. We 
create a map of the synchrotron spectral index in which 
the value at each pixel depends only on Galactic latitude 
(first 4 bins from Figure [3]). We then compare the fore- 
ground model derived using this spectral index map to 
a model derived using a uniform spectral index equal to 
the mean of the best-fit spectral index map. Allowing 
the spectral index to vary with latitude near the plane 
improves the x^ by 88 using only 4 additional parame- 
ters. 

We also test the assumption that the polarized dust 
distribution may be traced by the polarization angle de- 
rived from observations of polarized starlight. We repeat 
the fit in Eq. [51 replacing the dust polarization angle 7^ 
with the synchrotron polarization angle 7^ derived from 
the WMAP K-band data. The fit is marginally better 
(Ax^ = 9) when the polarized dust emission follows the 
magnetic field inferred from polarized starlight instead 
of the magnetic field derived from the synchrotron com- 
ponent. 

We may also use the polarization data to constrain 
emission from any additional component such as spinning 
dust. Direct comparison of the polarized maps yields 
spectral index (3s ~ —3.05 near the Galactic plane, con- 
sistent with synchrotron emiss ion but inconsistent wit h 
the inverted spectrum found bv lBoughn fc Pob"ei1 (|2006D . 
We test for polarized emission traced by the dust mor- 
phology by modifying Eq. [5] to include additional terms 

Qspm{n, = fspin{n)Td{n)Dspin{v) cos(27d(n)) 

Uspin{n, y) = fspin{n)Td{n)Dspin{v) ain{2jd{n)) (10) 

where /spin is the fractional polarization of the spinning 
dust component with f requency depend ence Dgpinii') 
peaking near 20 GHz ( Draine fc Lazarianlll998h . We re- 
peat the fit using either a spatially invariant synchrotron 
spectral index or the spatially variable map of polarized 
synchrotron emission derived from the K-Ka band com- 
parison (Fig. [21). We fit the spinning dust term within 
the same coarse pixels as the thermal dust emission, but 
make no assumption concerning the relative amplitudes 
of the thermal vs spinning dust components. In all cases, 
the fitted values are consistent with no spinning dust con- 
tribution to the polarization data, /spin = 0. Spinning 
dust contributes less than 1 percent of the observed po- 
la rization sig n al var iance in any band. 

iPage et al.l (|2007[ ) analyzed the WMAP polarization 
data outside the P06 mask to detect a cosmic signal at 
amplitude £{£+1)0^^/2^ = 0.086±0.029 ^K^ averaged 
over multipole moments ^ = 2 to 6. Their analysis used 
a foreground model derived by fitting template maps to 
the high-latitude portion of the WMAP Stokes Q and U 




5 10 15 20 25 30 35 
Multipole Moment t 



Fig. 8. — CMB power spectra outside the P06 mask, derived 
from Q and V band data after subtracting two different foreground 
models, (top) EE power spectra, (bottom) BB power spectra. 
The CMB spectra are robust with respect to different foreground 
models. 

data. This paper extends the foreground model to in- 
clude spatial variation and frequency dependence of the 
synchrotron spectral index along with spatial variation 
in the dust fractional polarization. The extended fore- 
ground model allows an improved understanding of the 
astrophysical foregrounds but does not significantly al- 
ter the cosmological conclusions. Fitting the variance 
of the residual maps as a superposition of instrument 
noise plus a cosmological signal yields a best-fit value 
= 0.38 ± 0.21 ^K^ for the CMB variance in thermo- 
dynamic temperature units, compare d to the expected 
value = 0.14 /iK^ derived from the iPage et al.l (|2007[ ) 
detection. 

The cross variance analysis provides a simple, model- 
independent limit to residual signals in the cleaned maps. 
Errors in the instrument noise model affect the quoted 
uncertainty but do not bias the fitted value. However, 
the variance analysis does not account for spatial corre- 
lations in the instrument noise. We have repeated the 
likelihood analysis by correcting the Q and V band po- 
larization maps using the pixel-based foreground model 
instead of the template model. Figure [51 shows the power 
spectra of the CMB signal outside the P06 mask, derived 
from the weighted sum of the cleaned cross spectra at Q 
and V ba nds on angu l ar sc ales up to the pixel size (c/ 
Fig 21 of IPage et all (|2007[ )1. The two cleaning tech- 
niques produce nearly identical results; the only change 
larger than the noise uncertainty is for the BB power 
spectrum at multipole moment £ ~ 2. We quantify the 
effect of different foreground cleaning on cosmological pa- 
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rameter estimation by jointly fitting the EE and TE data 
using the stand-alone likelihood code without fore ground 
marginalization as described in iPage et aTl (|2007( ) Table 
9. Changing the foreground model moves the fitted value 
for the optical depth from r = 0.092to;o3o to t = 0.098, 
a shift of roughly a fifth of the uncertainty. Systematic 
uncertainties in cosmological parameters caused by the 
foreground polarization are small compared to the noise 
levels of the WMAP 3-year data. 

5. CONCLUSIONS 

WMAP detects polarized emission over the full mi- 
crowave sky. The frequency dependence of this emis- 
sion is consistent with a superposition of polarized syn- 
chrotron and thermal dust emission. A pixel-by-pixel fit 
to the Stokes Q and U maps from the band-averaged 
3-year data separates the emission into synchrotron and 
dust components. Synchrotron emission dominates the 
polarization below frequencies of about 60 GHz. The 
spectral index of the synchrotron component, derived 
solely from polarization data, is -3.2 averaged over the 
full sky, with a modestly flatter index on the Galactic 
plane. Comparison of the polarized synchrotron emission 
to a maximum-entropy model of the unpolarized emis- 



sion shows typical fractional polarization of 2-4% near 
the Galactic plane, rising to 20% at latitude |6| > 50°. 
Prominent structures such as the North Galactic Spur 
appear as localized enhancements in the polarization 
fraction. We detect polarized dust emission with frac- 
tional polarization increasing from a minimum of 1% 
near the Galactic center to 6% at the anti-center. Dif- 
fuse emission from high-latitude dust is also polarized 
with mean fractional polarization 0.036 ± 0.011. This 
simple two-component model of the interstellar medium 
accounts for at least 97% of the polarized emission in the 
WMAP maps of the microwave sky. 
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